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C/5 ' Abstract. - We present an investigation of the electromagnetic instabilities which are trig- 

^ , gered when an ultra relativistic electron beam passes through a plasma. The linear growth rate 

O ■ is computed for every direction of propagation of the unstable modes, and temperatures are 

• ' modelled using simple waterbag distribution functions. The ultra relativistic unstable spec- 

rj , trum is located around a very narrow band centered on a critical angle which value is given 

'T^ • analytically. The growth rate of modes propagating in this direction decreases like k~^'^. 

The interaction of a relativistic electron beam with a plasma is a subject of relevance 
from many fields of physics ranging from incrtial confinement fusion [1] to some astrophysical 
K*" ' scenarios [2-4]. The linear analysis of the interaction reveals an highly unstable situation 
which has been investigated for a long time. Modes propagating along the beam are unstable 
within a certain range of wave-vector and form the electrostatic two-stream instability. The 
^p , so-called filamentation instability is found for modes propagating transversely to the beam 
f^ ' and finally, it can be proved that some modes propagating at arbitrary angle to the beam 
^^ I are also unstable [5]. As far as the growth rate is concerned, we see that it is eventually 
^^ ' a function of the parallel and perpendicular wave vector components. As long as the beam 
Y^ ', is not relativistic, the largest growth rate are found for wave vectors parallel to the beam 
which means that the two-stream instability dominates the linear evolution of the system in 



in 

o 



o 



^ I this regime [6] . The situation evolves when considering a relativistic electron beam. Because 
^H ^' relativistic electrons are harder to move in the direction of their motion than in the transverse 
O , direction, the two-stream growth rate is much more reduced than the growth rate of the modes 
propagating transversely, or even obliquely. If we denote 7f, the beam relativistic factor, the 
maximum two-stream growth rate is scaled like 7fj , the filamentation growth rate like 7^ 

rS I whereas the growth rate of the most unstable oblique wave vector is scaled like 7^ ' [5] . This 
j^ ' shows that oblique instabilities should dominate all the more than the beam is relativistic. 
The ultra relativistic limit is relevant for astrophysical settings such as high energy cosmic 
rays or gamma ray bursts production scenarios, for which some authors consider relativistic 
factors up to 10^ [2] and even 10^ [3]. These unstable oblique modes were first investigated in 
the cold limit (fluid approximation), and a temperature dependant treatment has only been 
developed recently [7]. As we shall see in this letter, accounting for temperatures in the beam 
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Fig. 1 - Numerical evaluation of the growth rate for a cold beam and a cold plasma, in terms of 
Z = kVij/ojpp. The growth rate is normalized to its maximum value 5o given by eqs. Q- Parameters 
are rib/np = 0.1 and 7b = 200. 



and the plasma results in a very narrow oblique unstable spectrum in the ultra relativistic 
limit. 

We thus consider a relativistic electron beam of velocity Vf, — Vtez, gamma factor jt, = 
1/(1 — V^/c^y^'^ and density nt passing through a plasma at electronic density Up ^ rib. Ions 
are supposed to form a fixed neutralizing background and we account for a plasma return 
neutralizing current [8] at velocity Vp such as n^Vf, = ripVp. The system is therefore charge and 
current neutralized, and we study the stability of harmonic perturbations ex exp(zk • r — ioj). 
We implement a 3D formalism using the linearized relativistic Vlasov equation as well as 
Maxwell's ones. Given the symmetries of the problem, the wave vector can be chosen within 
the {x, z) plan, and the dispersion equation reads [9] 
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in terms of the dielectric tensor elements 
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where /op and /of, are the equilibrium distribution functions for the plasma and the beam, m 
the electron mass, Wpp^b the electronic plasma frequencies for the plasma and the beam, and 
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7 = (1 + p^/to^c^)^/^. The dispersion equation (Q bridges between the two-stream modes 
at kx = and the filamentation ones at A:^ = 0. Because the former modes are longitudinal 
(k II E) while the later are purely transverse (k • E = 0), unstable oblique modes are neither 
longitudinal nor transverse and a fully electromagnetic dispersion equation such as 1^ is 
required to reach them. 

If we start considering a cold mono energetic beam entering a cold plasma, the resulting 
growth rate map is the one pictured on figure ^ We plot the growth rate normalized to the 
maximum oblique one in terms of the reduced wave vector Z = \<iVb/u)pp. In this cold limit, 
the maximum two-stream, filamentation, and oblique growth rates read 



One can check that the two-stream instability is negligible compared to the oblique and the 
normal modes. We also note on fig. ^ that the most unstable modes form a continuum at 
Zz ^ 0.9, starting from Z^ > 3. 

The introduction of temperatures modify the picture in various ways. As long as they 
remain "small" , they can be modelled through some simple waterbag distributions and the 
tensor elements Q can be calculated exactly [7]. It turns out that the resulting dispersion 
function has various singularities corresponding to various kind of preferential coupling be- 
tween a mode and some electron population. As the wave vector is more and more oblique, 
these singularities evolve and two of them overlap for a critical inclination angle 6c with the 
beam, which ultra relativistic expression is simply 

, n l + nb/Up 

tan0c = -i^7^, (4) 

vtpl Vb 

where Vtp is the plasma thermal velocity. At this junction, let us comment what is meant by 
"small temperatures" . As far as the plasma temperature is concerned, it simply means that 
the thermal velocity must remains small compared to the beam one. The present analysis 
is therefore valid for non-relativistic plasma temperatures, allowing for a 10 keV plasma for 
example. As far as beam temperature is concerned, transverse and parallel temperatures can 
be treated differently. Let us assume a waterbag distribution function defining a momentum 
spread around Pb — jb'mVb, 

f°b = ,„"'„ [e(p^ -Pb + Ptbw) - e(pz + Pb + Ptbw)] [e(p. + Ptb±) - e(p, - Ptb±)] , 

(5) 
where O is the step function. As can be seen on eqs. (0), the tensor elements are mostly de- 
pendent of the velocity distribution through the quantities p/7, and this is perfectly expected 
since the stability of a mode is a matter a wave particle interaction. Let us then evaluate 
the velocity spread {Av\\,Av±) corresponding to the momentum spread {Ptb\\,Ptb±) defined 
above. In the present ultra relativistic regime, we find 

Av« -- -^ , Av± r^ . (6) 

We recover for Av± the velocity spread corresponding to the momentum spread Ptb±- But 
A-yii is reduced by a factor 7^ so that the parallel beam temperature can be neglected in the 
very large 7b limit. 
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Fig. 2 - Numerical evaluation of the growth rate for a hot beam and a hot plasma, in terms of 
Z = kVt/ojpp. The growth rate is normalized to its maximum value 5o given by eqs. Q. Parameters 
are rib/np = 0.1 and 75 = 200. The parallel beam velocity has been set to 0, and every other thermal 
velocities are c/10. 



Considering a relativistic beam with 7b — 200, nb = n^/lO and neglecting its parallel 
thermal spread, we plot on fig. [21 the growth rate in terms of Z = kV^/wpp. The continuum 
of the most unstable modes observed on fig. ^ has turned into a very thin oblique region 
where the growth rate is higher than 90% of its cold oblique value 5o given by eq. Q up to 
Z ~ 50. Let us emphasize that unlike fig. Q] where the largest reduced wave vector plotted 
is Zz = 1.5 and Z^ = 5, we here extend the plot up to Z^ = 5 and Z^ = 50. Modes located 
above the critical angle are unstable up to Z^ ^ Jb/iVtp/Vb) + Z^ tan 0c [7]. We also note 
that the transition towards stable modes bellow the critical angle is very sharp so that the 
instability domain is almost rigourously bounded by the directions 9c and 7r/2. The critical 
angle obtained from eq. Q is arctan(ll) ^^ 84.8°, which corresponds very precisely to what is 
observed on the graph. In order to study the growth rate at large Z in the critical direction, 
we plot on figure O the maximum growth rate at angle 9c for Z up to 10"*. We observe a 
decrease with Z which can presently be fitted by 6/So — 2.1Z^^/^. It is interesting to notice 
that the maximum filamentation growth rate is here Sp — OAlSo so that the growth rate in 
the critical direction remains larger than the filamentation one until Z ~ 150. 

In conclusion, we have evaluated the growth rate of the electromagnetic instabilities for an 
ultra-relativistic beam passing through a non-relativistic plasma for any angle of propagation 
of the unstable modes. The parallel temperature of the beam can be neglected, and the 
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Fig. 3 - Numerical evaluation of the growth rate in the critical direction up to Z = 10^. Same 
parameters as fig. |21 The black points are the result of numerical calculation and the line represents 
5/5o =2.\Z-^I'\ 

unstable waves are found in a narrow domain comprised between 7r/2 and the critical angle Be 
which does not depends on the beam energy (see eq. 2)). Within this region, the most unstable 
modes are concentrated on a very narrow band extending around the critical direction, with 
a growth rate decreasing like Z~^/'^ oc k^^l'^ . 
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